Operator solutions for fractional Fokker-Planck equations 
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I. INTRODUCTION 

Ordinary derivatives account for the variation of a 
given function with respect to a given variable. Frac- 
tional derivatives have a more subtle meaning. We use 
throughout the Euler's definition of the fractional deriva- 
tive according to which the derivative of order a (0 < a < 
1) of a constant is indeed not zero, but 9" 1 = T{i-a) 
Their role in modelling physical phenomena is not 
intuitive and the treatment of the associated fractional 
differential equations (FDE) requires extreme care, not 
only from the mathematical point of view. The gener- 
alization of a relaxation equation, with a constant force 
term, to a fractional form reads [1, Q 



r(l - a) 



P + f, 



(1) 



where P = Pa{0) is the initial condition and k is a con- 
stant. Eq. ([1]) is an ath order FDE with / being the 
non-homogeneous part. The term with P is not a gen- 
uine non-homogeneous contribution, but it accounts for 
the nonvanishing of a constant under fractional deriva- 
tive. We use in Eq. ([T]) the Euler definition of fractional 
derivative because it appears most suitable to treat the 
dynamical behavior governed by the FFP equation we 
will discuss later. The problem ([T]) is mathematically 
well defined. The apparent singularity at t = can be 
removed by multiplying both sides of the equation by 
dl~", thus getting 



dtPc.it) = dl-^{-nP^{t) + f). 



(2) 
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The notion of stationary solution is well defined for 
an ordinary relaxation differential equation, but not for 
its fractional counterpart. In common terms stationary 
means that the solution is no more sensitive to time vari- 
ations and, hence, its (ordinary) time derivative is zero. 
The notion should be revised for FDE, in accordance with 
the order of the derivative. The solution of Eq. ((TJ reads 



(3) 



where Ea,fj{z) ~ 'Yl'^^o^'' l^i'^''' P) '^^ ^^"^ modified 
Mittag-Leffler function and reduces to its ordinary case 
for /5 = 1, (z) Q. 

The solutions, plotted in Fig. [T] for different values 
of a, do not display any long time stationary behavior. 
Quasi-stationary behavior is reached for a approaching 
the unity. For large t we find Pa{t) oc t"^^, for which the 
ath derivative is vanishing. We can therefore conclude 
that the solution reaches a- derivative stationary form. 



II. FRACTIONAL FOKKER-PLANCK 
EQUATIONS AND EVOLUTION OPERATORS 

The Eq. (2) can be generalized to the following partial 
differential equation 



d{'Fo,{x,t)^LFpF^{x,t) + 



t~ 



r(l - a) 



7(x), (4) 



which has been shown to be tailor suited for study of 
problems of anomalous diffusion The initial condi- 
tion is Fa{x,0) ~ j{x). From mathematical point of 
view Eq. (|4]) is a well-posed Cauchy problem and it is 
the two- variables generalization of Eq. ([Ij. In Ref. Q 
Eq. (U) has been used to model the continuous time ran- 
dom walk with the inclusion of effects of space dependent 
jump probabilities and Lpp denotes the Fokker-Planck 
(FP) operator involving the spatial derivative dx- The 
presence of the term with 7(0;) ensures that Eq. (U) is 
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FIG. 1: Logarithmic plot of Pa{t), (see Eq. ((Sj), for k = 1, 
P = 1, / = 2 and a = 1/4 (I), 1/2 (II), and 5/6 (III). 



well defined and describes a process preserving the norm 
of the distribution Fa{x,t)^ when the time evolves. For 
any function h(x) its average value over Fa{x,t) is de- 
fined as 



{h{t))c, = / h{x) Fo,{x,t)dx . 



(5) 



The formal solution of Eq. is obtained by using 
an extension of the evolution operator formalism, intro- 
duced by Schrodinger, therefore getting 

F^{x,t) = Uc.{t)j{x), Ua{t) = E^{t"LFp). (6) 

Below we shall apply Eqs. ^ to three different ver- 
sions of Fokker-Planck operators Lpp. Limiting for the 
moment the discussion to L pp — kd^, where k is the 
generalized diffusion constant, Eq. (j4]) can be interpreted 
as the fractional version of the heat equation @ and its 
solution reads 



{kt^y 



F^ix,t) =J2wT_^ , 
^ Ffl -far) 

r—O ^ ' 

which for 7(x) — x" (n e Z), gives 



[9^7^], (7) 



V2] 



(n-2r)!F(l + ar)' 
(8) 

which are polynomials in x. For a = 1 they are 
known as heat polynomials @. Any initial function 
li^) = X^^o'-"-^" allows therefore a solution of the 
fractional heat equation as the following expansion 



OO 

E 

n=0 



(9) 



As in the case of conventional heat equation, the se- 

(2) 

ries in terms of fractional heat polynomials aHn are of 
limited usefulness since it converges for short times only. 
As an example, for a = 1 and 



1 



' 27r (Jx 



(10) 



the convergence is limited to t < (j^/(4 /c). The use of the 
Gauss- Weierstrass transform Q provides solutions with a 
well behaved long time behavior and therefore we look for 
an analogous transform relevant for the fractional case. 

We make therefore the assumption that such a trans- 
form exists and that we can write 



na{s,t) e"" ds 



(11) 



with na{s,t) being not yet specified functions. The evo- 
lution operator in Eq. ^ can therefore be written as 



Uc.it) = / 7i„(s,0C/i(s)ds, J7i(t) 



„tLB 



and, equivalently, 

Fo.{x,t) 



na{s, t) Fi{x, s) ds 



(12) 



(13) 



holds. Eq. provides the link between fractional and 
ordinary Fokker-Planck equations through the knowledge 
of Ua {s,t). This equation, specified to the case of Eq. ([8]) , 
leads to the following relation 



nc.(s,t)iHi^\x,s)ds, (14) 



which yields, direct consequence of Eq. PT]) 



n'ais,t) — - ds 
ml 



r(l + a m) 



(15) 



According to Eq. (16) in [ll|, the functions na{s,t) can 
be identified as 



rials, t) 



t 



1 



5o 



t 



(16) 



The functions ga{z) are the one-sided Levy stable dis- 
tributions, recently obtained for a rational in d, [H. For 
related considerations compare Il0| . The Eq. (jl2p is simi- 
lar to the one reported in Refs. [5| and ll|. The meaning 
of the fractional evolution operator is that the solution of 
the fractional Fokker-Planck (FFP) equation of order a 
is known whenever that the ordinary case, a = 1, is avail- 
able. By simple manipulation of the previous equations 
(see Eqs. pT|) and we can also conclude that 



ilpit) ^ / np/a{s,t)iJa{s)ds, I3<a. (17) 
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Therefore, the solution of the FFP equation of order 
P can be derived from its a counterpart by a self- 
reproducing procedure. It should also be noted that, for 
a ^ 1, Ua{t2 + h) ^ Ua{t2)U{ti). The evolution at 
different times ^2 > ii is therefore given by 

/•OO 

Uc.{t2)Ua.{ti) = / n„(si,ti)dsi (18) 

POO 

X / na{s2,t2)ijl{si+ S2)AS2. 

Jo 

This formula turns out extremely useful to deal with suc- 
cessive approximations, when the nature of the FP op- 
erator does not provide any close form for the operator 
Ui{t). The functions na{x, t) defined in Eq. (|16|) turn out 
to be, for a = fc = 2, 3, . . ., solutions of general heat 
equations dt ui/i^{x, t) — (—1)'"" Ui/k{x, t) with the ini- 
tial condition 0) = 5(x). These heat equations 
have been also obtained in jl2l | from purely probabilistic 
arguments. The case of Ua{x,t) for a = l/k and I > 1 
will be the subject of a forthcoming study. 



Glaisher form a shift term in the Gaussian provided by 
Tt/{mor]a). The solution for different values of a and 
t = 2 are given in Fig. [21 and the first moment of the 
distribution is, see Refs. Q and [T]| : 



FJx, 2) 




1 2 3 4 5 6 7 

X 



III. SPECIFIC EXAMPLES 

We can now apply the wealth of the operator tech- 
niques known for the conventional FP equation, to solve 
its fractional version. For instance, starting with Gaus- 
sian initial condition of Eg. (j lO p, we evaluate Ui{s) ^{x) 
with the Glaisher formula Il3l. 11411 and obtain 



Uiishix) 



1 



\/27r (Ta 

X exp 



2kis 



X 

'2ai 



-1/2 



2kis 



(19) 



which, according to formula ([5]), gives (a;^(s))i = 
cr^ (l + and, by using Eq. P^ . allows us to con- 

clude that the a- and ^-dependent variance of x is given 
by 



poo 

aljt) = / n^{s,t){x^{s)) 
Jo 



id.s 



(20) 



2kt° 



" r(l-fa)' 

Note that for j{x) = 5{x), we have formally cr^ = and 
Eq. (|20|) reproduces the defining equation of subdiffusive 
behavior. 

We now move on to more general Fokkcr-Planck 
operators. We start by considering the operator 
Lpp = kd^ + [T/iruQija)] dx, where the second term 
is due to the action of a constant force rja is fractional 
friction coefficient, and mg is the particle mass. The so- 
lution of our problem can be written by adding to the 



FIG. 2: Plot of the solution of Eq. @, Fa(x,t), with FP oper- 
ator Lfp = dx + 2da: and initial condition 7(2;) = ^"^ , 
for t = 2 and Q = 1/4 (I), 1/2 (II), and 3/4 (III). 

The operator Li?p — [al d"^ + dxx) is used in storage 
ring physics to model the effect of diffusion and damp- 
ing (r is the damping time) of the electron beam due 
to the synchrotron radiation emitted by the electron in 
the bending magnets of the ring [l5| . tXe is the variance 
of so-called equilibrium distribution (see below). The 
two processes, (damping and diffusion), yield eventually 
a stationary solution in the conventional case. Such a 
condition does not exist for the fractional version. 

The evolution operator Ui{t) associated with the last 
FP operator can be written in a simple form. By setting 

indeed A = and B = 2-dxX we obtain A, B = 

^a^A, so that the use of conventional operator ordering 
methods yields [ll] 



Uiit) 



exp 



1 - e-^*/" 



A 



(22) 



In the case in which the initial distribution is the Gaus- 
sian the fractional evolution will be characterized by the 
following variance 



with 



{x\t)), 



na{s,t) {x'^{s))ids 



(23) 



(24) 



where a is the variance of the initial distribution. In 
Eq. ((221) {x'^{t))a is obtained by replacing, according to 
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Eq. (HI]) in the second term e"''*/'^ with Ea{-'it°' /t). 
(Note that the physical dimension of the damping time 
T is [t"]). In Fig. [3] we reported the {x'^{t))a and, as ex- 
pected, equihbrium conditions in the conventional sense 
is not reached. The plot shows however the onset of anal- 
ogous regimes after the knee-shaped decrease. This is a 
consequence of the fact that for increasing time the sec- 
ond term in Eq. (|24p becomes dominating with respect 
to the first, and all solutions approach the a-derivative 
stationary form. 




FIG. 3: Double logarithmic plot of {x^{t)}c, (see Eq. ((23)) '). 
for (7 = 2, cte = 1, and a = 3/5 (I), 4/5 (II), and 9/10 (III). 



IV. DISCUSSION AND CONCLUSION 



In such a context the Levy stable distribution function is 
going to play a role in the theory of FFP of type Eq. ([25l) . 
The use of the evolution operator yields a formal solution 
of the type Ga{x, t) = e~^*^° h{x). The series expansion 
of the exponential may have a limited use only, we look 
therefore for a more useful representation of the evolution 
operator. The use of the identity [1] 



e-"^" = 5a(e) exp(~ai/>e)dC (26) 
Jo 

is the naturally suited choice, so that we find 



G^{x,t) = (At) 



-I/O 



ga 



(At)i/c 



h{a)A(T. (27) 



This technique (albeit limited to the case a = 1/2) has 
been applied to the study of the relativistic heat equation 
{dt Gi/2{x, t) = — -^/l — 92 Gi/2{x, t)) in [l^ and appears 
a very promising tool in further applications, possibly 
involving relativistic quantum mechanics. 

Finally, we emphasize that the solutions of Eq. ^ for 
1 < a < 2 can also be obtained with the help of the 
evolution operator and of two-sided Levy stable distri- 
butions obtained in 3. The form of Eq. ^ has to be 
however modified as then the 1^-^(3,1) function has to be 
replaced by its two- variable counterpart discussed in [9| . 
In this context we refer to Eqs. (5.21) and (5.22) of [l2[ 
where the two-sided case is studied for the heat-type FP 
equation. 

The different topics touched on in this paper have 
shown that the combined use of techniques from vari- 
ous fields (including statistical mechanics, theory of frac- 
tional calculus, ordinary quantum mechanics, etc) offers 
the appropriate tool to study new phenomena in the the- 
ory of anomalous diffusion. 



We can also consider the case of partial fractional dif- 
ferential equations in which the fractional derivatives is 
acting on the spatial coordinates. From the mathemati- 
cal point of view we have the following Cauchy problem 

dtGc,{x,t) = -Xd^Gc,{x,t) (25) 
Ga{x,0) = hix), 
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